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We investigate quantum information processing and manipulations in disordered systems of ultra- 
cold atoms and trapped ions. First, we demonstrate generation of entanglement and local realization 
of quantum gates in a quantum spin glass system. Entanglement in such systems attains significantly 
high values, after quenched averaging, and has a stable positive value for arbitrary times. Complex 
systems with long range interactions, such as ion chains or dipolar atomic gases, can be modeled 
by neural network Hamiltonians. In such systems, we find the characteristic time of persistence of 
quenched averaged entanglement, and also find the time of its revival. 



Successful implementations of quantum information 
processing (QIP) in atomic, molecular, or solid state sys- 
tems typically demand very rigorous control of such sys- 
tems This concerns both few qubit systems such as 
the Cirac-Zoller computer with ions or photons 
as well as atomic gases in optical lattices |4|. Despite 
a lot of progress, the demanded control in such systems 
is nowadays very hard to achieve Recently QIP in 
systems with a limited knowledge of the parameters has 
also been proposed Q. 

At the first sight, what we propose here sounds like 
contradictio in adjecto: QIP in quenched disordered or 
complex, ergo hardly controllable, systems. However, as 
we have recently shown, one can create controlled dis- 
order in atomic gases in optical lattices and study, in 
an unconventional way, Anderson and Bose glasses in a 
Bose gas |2| , or spin glasses with short range interactions 
in Fermi-Bose, or Bose-Bose mixtures ||. Using linear 
chains of trapped ions 0, or dipolar atomic gases [Io| . 
it is possible to realize complex spin systems with long- 
range interactions that may serve as model for classical 
and quantum neural networks [ll| . 

Disordered systems offer at least two possible advan- 
tages for QIP. First, they have typically a large number 
of different metastable (f ree) energy minima, as it hap- 
pens in spin glasses (SG) [l^. Such states might be used 
to store information distributed over the whole system, 
similarly to neural network (NN) models The in- 

formation is thus naturally stored in a redundant way, 
like in error correcting schemes ^ij. Second, in disor- 
dered systems with long range interactions, the stored 
information is robust: metastable states have quite large 
basins of attraction in the thermodynamical sense. 

We address here the simplest fundamental questions 
concerning QIP in disordered or complex systems: (i) 
Can one generate entanglement in such systems that 
would survive quenched averaging over long times? (ii) 
Can one realize quantum gates with reasonable fidelity? 



Here we answer both questions affirmatively considering 
both short and long range disordered systems. 

First, we consider a short range disorder Ising Hamil- 
tonian, the so-called Edwards-Anderson (E-A) model of 
spin glasses which can be straightforwardly implemented 
using atomic Bose-Fermi, or Bose-Bose mixtures in opti- 
cal lattices |il Il5|. We address the generation and evolu- 
tion of nearest neighbor (nn) entanglement in this model. 
In the short range Ising model without disorder, it is 
possible to create cluster and graph states (i.e. entangle- 
ment) starting from an appropriate initial product state 
(lif. Here we show that, while the disorder averaged den- 
sity matrix of two neighboring spins remains always sepa- 
rable, the disorder averaged entanglement (quantified by 
logarithmic negativity [rjj) converges with time to a fi- 
nite value. The generation of entanglement 0] as well as 
its evolution for arbitrary times in an Ising model without 
disorder but with long-range interactions, has also been 
addressed in Rcf. There it was suggested the possi- 
bility of applying similar ideas to disordered systems. We 
show also that the quantum single-qubit Hadamard gate, 
can be realized in such system with significant (disorder 
averaged) fidelity. 

Secondly, we consider complex systems with long range 
(1/r 3 , or 1/r 2 ) interactions, that can be realized for in- 
stance, in linear ion traps, using either local magnetic 
fields, as proposed by Wunderlich and coworkers or 
by appropriately designed laser excitations . The cor- 
responding Hamiltonian can be mapped into an I sing 
Neural Network (NN) model with weighted patterns [13j . 
Those patterns can be used as qubit systems, with the 
information distributed over the chain. One can also in- 
clude external parallel, or transverse fields in the model. 
We show that in such system, it is possible to generate 
long range bipartite entanglement that undergoes a series 
of collapses and revivals |2flj , whose times are found an- 
alytically. Finally we study also bipartite and tripartite 
entanglement dynamics in an infinite range Ising model 
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without disorder. 

Let us start with the Edwards- Anderson spin glass 
model described by 
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Here cr| denotes the Pauli operator at the kth. site, and 
Jy's describe nn couplings for an arbitrary lattice. In 
the E-A model these couplings are given by independent 
Gaussian variables with mean J and variance a 2 . Start- 
ing from a pure product state of the form \^f) = ]J i 
where |±> = (|0>±|l»/ v / 2 H3 , we evaluate the entangle- 
ment after a finite time, where the density matrix is given 
by p(t,{Jij}) = exp{-iH E -At}\^){^\exp{+iH E -At}. 
The reduced density matrix for a nn pair is obtained by 
tracing over all other sites. For instance, the reduced 
density matrix for a 2D square lattice is given by 
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:- iJl2t / 2 |cOs(J 13 V2cOs(Ji5t/2)cOs(Ji7V 2 )|01>(H| 



i(J24t/2)cos(J 26 t/2)cos(J28t/2)|10)(ll|} 
cos( Ji3t/2) cos( J\^t/2) cos(Jnt/2) cos( J 2 4i/2) 
i( J 26 t/2) cos( J 28 i/2) { |00) <11| + |01) (10| } + 



h.c. 



where 1 is the identity operator and the indices 3 ... 8 
enumerate the six neighbors of 1 and 2. A similar ex- 
pression can be obtained for the ID lattice. In both 
cases, the averaging of the reduced state over Jy 's (equiv- 
alent to reducing the average Q\i{t, { Jy})) is separable. 
Note, however, that as always in physics of disordered 
systems, if we are interested in typical values of physical 
quantities such as free energy, entanglement, etc., we are 
obliged to perform a "quenched" average, i.e. first cal- 
culate t he q uantity of interest and then average ^3] ( see 

also mm!)- 

To study entanglement, we use the logarithmic neg- 
ativity (LN) ^?|]- The LN of a bipartite state pab is 
defined as E LN {p AB ) = log 2 ||p^||i, where ||.||i is the 
trace norm, and p^ denotes the partial transpose of pab 
with respect to the ^4-part [23] ■ Note that pij (t) acts on 
C 2 Cg>C 2 . Consequently, a positive value of the LN implies 
that the state is entangled and distillable [23I |3| , while 
Eln — implies separability [23|. 

The entanglement in the spin glass model turns out 
to be an even function of the couplings. The temporal 
behavior of EL^(t) in a 2D square lattice is shown in 
Fig. H for two different cases of disorder: with frustra- 
tion and without it. For J = 0, a 2 = 1, the system has 
randomly ferro- (J > 0) and antiferro-magnetic (J < 0) 



interactions and is strongly frustrated; E^Mit) is rapidly 
damped to a constant, and does not show any oscilla- 
tions. This behaviour differs from the non-frustrated 
case J = 5, a 2 = 1, when EL^it) exhibits oscillations 
with frequencies ~ 1/J. For short range interactions, 
the next-nearest neighbor entanglement vanishes, even 
before the averaging, for both ID and 2D. To understand 
why entanglement converges in time to the same finite 
value in both the frustated and non-frustated cases, no- 
tice that as long as the distributions Jy 's are sufficiently 
well-behaved, Jijt/2 corresponds to a uniform distribu- 
tion over [0, 2ir] for large enough t. 
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FIG. 1: Temporal behavior of nn averaged entanglement in a 
2D spin glass model, starting from IT |+}-. For a model with 
frustration (J = 0), Eln it) converges quickly to a constant 
value (red curve). For a non-frustrated case (J = +5), -ELJv(i) 
exhibits damped oscillations (blue curve), converging to the 
same value « 0.0154, as reached in the frustrated case. Stan- 
dard deviation for t — » 00 is ~ 0.0704. It is interesting to note 
that the dynamical behaviour of Eln depends on J, although 
at large times, they all converge to the same value. The same 
behavior is encountered in the ID case, even though there is 
no frustration in that case. 



We have calculated the nn entanglement for the fol- 
lowing lattice configurations: ID chain, 2D honey-comb 
lattice, 2D square, 3D cube, where any given pair of 
neighboring lattice sites has al = 2,4, 6, 10 neighbors re- 
spectively. For time large enough, our numerics reveal 
that bipartite entanglement decays exponentially with 
the number of neighbors. Such behaviour can be repro- 
duced analytically by considering the volume of the set 
of separable states (see e.g. |2aj). giving an upper bound 
on nn entanglement that depends exponentially on d. 
Some algebra shows that if the state p^ (t) is entangled, 
then £i=i$ < (3 - 4i? 2 )/2, where the & = Jyt/2's 
are state parameters varying from to 2-7r, and R is the 
radius of the separable ball in the c?-dimensional space. 



The volume of this hypersphere is Vd = Sd 



/d, 
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where Sd = 2tt^ /T(d/2). Due to the periodicity involved 
implicitly in Pij(t), there are 2 d — 1 such hyperspheres. 
Considering all states in this volume to have unit en- 
tanglement, the average entanglement at long times is 
£d = Vd(2 d — l)/(27r) d . As an example, for the case of the 
2D lattice (for which d = 6), at long times, the actual en- 
tanglement is w 0.0154, while £ 6 w 0.0221. Although the 
bipartite entanglement vanishes with increasing number 
of neighbors, one can expect the multipartite entangle- 
ment to be non vanishing due to the fact that the volume 
of separable states is "super-doubly-exponentially small" 
with increasing number of parties [2(| . 

We show now that spin glasses allows also to imple- 
ment quantum gates. We focus on the Hadamard gate, 
which transforms the computational basis into a comple- 
mentary basis: |0) — > |+) and |1) — > |— ). To implement 
the Hadamard gate, assume that the computation is per- 
formed in a spin lattice, and the particles 1 and 2 are a 
part of it. We assume that at a certain time, particle 1 
is in an arbitrary state a|0) + 6|1), where \a\ 2 + |6| 2 = 1, 
and we let system evolve according to the Hamiltonian 
He-a for a suitable duration of time, before perform- 
ing measurement on particle 1 (in a suitable basis). For 
J — 5, a 2 — 1, particle 2 attains the Hadamard rotated 
state a\+) +b\— ), with quenched averaged fidelity greater 
than 0.85. One can increase such fidelity by increasing 
the number of spins, and employing assisted measure- 
ments. Note, that if we try to prepare the Hadamard ro- 
tated state using the classical information obtained only 
from the measurement of particle 1, the fidelity is only 
2/3 [13. 

Let us now move to a long-range interactions spin 
Ising model, described by the Hamiltonian Hi r = 
W^ij Jij cr i cr j' wnere N is the total number of spins. 
Such models can be realized with trapped ions 1 1] , where 
J u = 12 tl Q£i/*l> wit h Q (A„) describing the phonon 
eigen-modes (eigen- frequencies) . Here we consider two 
extreme cases. First, we take Ai = l,£j = constant 
Vi, A M — > oo for fi > 2, so that the interactions are 
ordered, and the Hamiltonian is Hi ro = jjS 2 , where 
S = J2iLi a i- Secondly, we consider the case when 
A,, = 1 for all p, when the Hamiltonian becomes Hn n = 

jf E£=i EJUi & } d j) v?<?j- This is the Hopfield model 
of a neural network with Hebbian couplings 0] . Here p 

is the number of "patterns" of the neural network, and 

(i) 

the patterns are described by random variables = ±1, 
each with probability h. As in the case of short-range 
interactions, we take the initial state of the evolution 
as \^f) = n^ x |+)j, and study the dynamics of entan- 
glement for ordered and disordered Hamiltonians. We 
provide an efficient method to analytically compute the 
evolved state of any number of patterns and any number 
of spins. 

Consider first the case of the Hamiltonian Hi ro . 
We can write the evolution operator exp(—iS 2 t/N) as 



/ d(jjexp((i/N)uj 2 + S\/i(—2i/N)uj), up to aconstant fac- 
tor. Applying now this unitary to the initial state |^), we 
find any two-party state Qi^i^) = trfc/1,2 p lro {t) of such 
system and compute the entanglement quantified by the 
LN. (This method can be also applied to find multipar- 
tite evolved states). In Fig [21 we plot the entanglement 
(as quantified by LN) of £>i2 (i) ; with respect to time, 
as well as N. The figure shows revivals of bipartite en- 
tanglement, that occur on the time scale tr ~ N, and 
persist on the time scale tq ~ V~N (collapse time). 




FIG. 2: Generation of entanglement of bipartite states £?i2°(i) 
with respect to time and number of spins. Collapses and 
revivals of the entanglement are clearly depicted. 



As depicted in Fig. |3 there are large ranges of 
time, for which the bipartite state is separable. In- 
terestingly, this range of separability can be reduced, 
considering entanglement of the tripartite evolved state 
p l {2 3 (t) in a bipartite cut. Although the interactions in 
Hi ro are long-range, they are ordered, so that p&tf) 
and Pi23(i) takes a relatively simple form. Amazingly, 
the same method applies for Hnn, where the interac- 
tions are both long-range and disordered. Despite its in- 
creased complexity, we can still use the technique for the 
evolution operator exp(— iH^Nt), that was used in the 
case of Hi ro . Specifically, we replace in exp(— iH^^t), 
the operator exp(—iS 2 t/N) by f duj^exp^i/N)^ 2 + 

S^Vt(—2i/N)Uf 1 ), for every p, where 5 M = Eilid ' of. 
Applying this operator to our initial state, we find that 
the A^-particle state at time t is 



where r„ 



|-{e- 2< ^5: M €?V ( ./^(|o)(i|) i +h.c.}] J 

., = + a,/, = w M - uif/, with (jl = 1, 



(3) 
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After tracing out all except particles 1 and 2 we obtain: 

e ™(t) = i/4{(|oo)(oo| + |oi)(oi| + |io)(io| + |ii)(ii|) 

+ [e" 4it £ , e'4 2) /iv 2 CQa{4t J- $ffi/N)\00) (01| 

+ n^ li2 cos(4i tfkP/N^OO) (lof 

+ n^cosC^^ 1 ) +4 2 ))/AT)|00)<ll| 

+ n^ 1)2 cos(4^«(#> -e«)/iV)|oi)(io| 

+ [ e <* (n^ 1)2 cos(4< £ 101) (11| 

+ n i#1 , 2 cos(4t^^ ) ^ 2) /iV)|10)(ll|)] +h.c.]}. (4) 



For N large, and 
expression can be 
fact 



t/N small, 
simplified 



that n^ 1)2 cos(4t £ M sjPtP/N) 
E^i, 2 log e |cos(4iE Al ^ ) d 2) /^)l 



the 
using 



above 
the 



exp 

exp[- 

r i - 
•V 
fore, 



-(8i 2 /^ 2 )E^i, 2 (E M 4) 2 ] ; where for all i, 
- +1 or —1 with probability 1/2 each. There- 
for large N and small t/N, we have that 
n^i,2 cos(4f^^£^£^/iV) self-averages to the value 
cxp[— (8t 2 p/N)], so that after time t ~ y/N/p, all 
the off-diagonal elements of the state Q^ 2 N (t) become 
vanishingly small. Therefore, for the first time, nearest 
neighbor entanglement in the evolved state appears and 
persists for times of order tq ~ y/N/p. However, there 
are repeated revivals in entanglement, with the period 
being tr ~ ttN/2 for odd p, and tr — nN for even p. 
Note, that the period of revivals is independent of the 
number of patterns in the model (cf. [22j|). 

Summarizing, we have studied disordered and complex 
spin systems with short-range and long range interac- 
tions that can be realized with trapped atoms or ions. 
We have shown that in both cases it is possible to gen- 
erate quenched averaged entanglement over long times. 
In the case of short range interactions, we considered 
Edwards- Anderson model in ID and 2D square lattice. 
We have shown that in such disordered system, it is pos- 
sible to implement also distinctly quantum single-qubit 
gates with high fidelity. We have also demonstrated that 
it is possible to generate entanglement in the spin system 
with long range interactions, corresponding to the Hop- 
field neural network model. We have shown that in such 
case, entanglement exhibits a sequence of collapses and 
revivals. 
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